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O ^ Abstract. From the inversion of a time series of high resolution slit spectrograms obtained from the quiet sun, the spatial and 

temporal distribution of the thermodynamical quantities and the vertical flow velocity is derived as a function of optical depth 
(logT) and geometrical height (z). Spatial coherence and phase shift analyses between temperature and vertical velocity depict 
the height variation of these physical quantities for structures of different size. An average granular cell model is presented, 
showing the granule-intergranular lane stratification of temperature, vertical velocity, gas pressure and density as a function 
of log x and z. Studies of a specific small and a specific large granular cell complement these results. A strong decay of the 
temperature fluctuations with increasing height together with a less efficient penetration of smaller cells is revealed. The T - 
T coherence at all granular scales is broken already at log T = — 1 or z~ 170 km. At the layers beyond, an inversion of the 
temperature contrast at granular scales > l"5 is revealed, both in log T and z- At deeper layers the temperature sensitivity of the 
, H~ opacity leeds to much smaller temperature fluctuations at equal log T than at equal z, in concordance with Stein & Nordlund 

^ — ' (1998). Vertical velocities are in phase throughout the photosphere and penetrate into the highest layers under study. Velocities 

at the largest granular scales (~ 4") are still found even at log T ~ —2.8 or z ~ 370 km. Again a less efficient height penetration 
of smaller cells concerning convective velocities is revealed, although still at logT~ —2 or z~280km structures > 1"4 are 
detected. A similar size distribution of velocity and temperature structures with height provides observational evidence for 
substantial overshoot into the photosphere. At deep photospheric layers, the behaviour of the vertical velocities reflected in 
simulations is for the first time qualitatively reproduced by observations: Intergranular velocities are larger than the granular 
ones and, both reach extrema, where the granular one is shifted towards higher layers. 
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1. Introduction 1990, 1991a, 1991b, 1991c; Hanslmeier et al. 1991a, 1991b, 

1993, 1994, and references therein), but without filtering prop- 

During the last two decades the vertical structure of the solar erly for oscillationSj which strongly mask the convect i ve sign a- 

photosphere has been investigated by different methods. ture 

Statistical analyses of a few parameters like continuum in- 
tensity, line core intensity, velocities (derived from Doppler Altrock et aL ( 1984 ) and Nesis et al - ( 1999 ' 2001 ) followed 
shift measurements in the line core), equivalent width and full the evolution of selected granules during short periods of time, 
width at half maximum are often used, in addition to line bi- Keil ( 1984 ) and Johanneson et al. (1992) obtained 2D spec- 
sectors, in order to characterise the shape of photospheric lines tral ima S es derived from 1D spectrograms by scanning a cer- 
and to infer important physical parameters. ^ area of the solar surface - Collados et al. (1996) obtained 

. . t . .. , , , , , t . , improved data using the Correlation Tracker developed at the 

Most investigations have been based on correlation and r b r 

. , , • . , • ,. j . ■• • w,r>i\ Instituto de Astrofisica de Canarias. 
spectral analysis techniques applied to one-dimensional (ID) 

slit spectrograms (Nesis et al. 1988, 1992, 1993; Komm et al. 2D spectrograms have been obtained using a Multichannel 

Subtractive Double Pass (MSDP) spectrograph (Roudier et al. 

Send offprint requests W. K. G. Puschmann 1991; Espagnet et al. 1993, 1995) or Fabry-Perot interferome- 
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ters (Salucci et al. 1994; Bendlin & Volkmer 1993; Krieg et al. 
2000; Hirzbergeret al. 2001). 

On the other hand, one of the main issues of solar atmo- 
spheric studies is to determine in a more accurate way the 
stratification of physical quantities throughout the solar pho- 
tosphere. This can be done by different techniques. 

One is to assume a certain model atmosphere and to ob- 
tain through forward modelling synthesised observable param- 
eters which can directly be compared with observations. See 
to the numerical simulations of solar convection by Steffen 
et al. (1989, 2003), Stein & Nordlund (1989), Gadun et al. 
(1997, 1999, 2000), Freytag et al. (2002), and Wedemeyer et 
al. (2004). 

Another technique is to invert the observations in order to 
obtain a model atmosphere through iteratively comparing syn- 
thetic observables with real ones. Inversion methods address 
the problem of the determination of the stratification of physi- 
cal quantities by simultaneously varying all (or at least as many 
as possible) relevant physical quantities in a parameterised at- 
mospheric model in order to obtain iteratively the best fit to the 
observed line profiles. 

The first semi-empiric 2D model of the solar atmosphere, 
including a physical model of granulation, was developed by 
Nelson (1978). Later Skumanich & Lites (1987) obtained a 
model using their inversion code based on the Milne-Eddington 
approximation. Several inversion methods have been devel- 
oped at the IAC. Ruiz Cobo & del Toro Iniesta (1992) pre- 
sented the inversion code SIR (Stokes inversion based on re- 
sponse function), dealing with arbitrary stratifications of phys- 
ical quantities in the photosphere. Although the SIR code was 
applied to deduce the structure of granules and intergranular 
lanes in the sun (Rodriguez Hidalgo et al. 1996, Ruiz Cobo et 
al. 1996), the results have never been published in cited jour- 
nals due to three fundamental criticisms: The data have not 
been filtered for oscillations, the spectral lines used for these 
studies have not been corrected for departures from LTE, and 
one of these lines shows strong magnetic sensitivity. A very 
similar code has been applied on low spatial resolution but high 
signal to noise ratio spectra (Frutiger et al. 2000). 

The state-of-the-art knowledge of the granular phe- 
nomenon can be described in the following way: Convective 
overshoot of the plasma from the solar convection zone into 
the photosphere forms a pattern of bright cellular elements 
showing upwards motion - the granules - surrounded by a 
network of dark intergranular lanes, where down-flow mo- 
tions are observed. According to correlation analyses (Deubner 
1988; Salucci et al. 1994; Espagnet et al. 1995) the horizon- 
tal temperature fluctuations associated with these motions de- 
crease rapidly with increasing height until they vanish. Only 
the largest granules (> l."5) contribute to the brightness pat- 
tern observed above where the brightness contrast is inverted. 
Note that the values of the height where the temperature fluctu- 
ations vanish varies a lot in the literature. One finds values from 
60 km reported by Kneer et al. (1980), 60- 90 km (Espagnet et 
al. 1995) 170 km (Hanslmeier et al. 1993; Kommet al. 1991a) 
up to 270 km where Bendlin & Volkmer (1993) detected bright- 
ness signatures of the granulation. The values may differ from 
each other due to differences in the method used to establish 



the geometrical height scale (transformation from T to z) and in 
different methods to filter the oscillations. The vertical velocity 
field persists in the upper layers, but waves (acoustic and grav- 
ity) and turbulent motions may also contribute to the velocity 
fluctuations (Deubner 1988; Deubner & Fleck 1989; Salucci et 
al. 1994; Espagnet et al. 1995). 

In many of the studies above, based on the analysis of line 
parameters, specific positions in the spectral lines have been 
directly associated with geometrical heights by means of the 
'formation heights'. This is not a reliable concept since in real- 
ity, the information from different atmospheric layers is mixed 
in the spectral lines. Therefore, it is more appropriate to study 
the response functions to a physical parameter, which show the 
variation of the emergent intensity at a given wavelength pro- 
duced by a differential change in the specific parameter at a 
given photospheric layer (Sanchez Almeida et al. 1996). 

In this series of papers a complete study of the vertical 
structure of the solar photosphere will be performed. After a 
first attempt with a global correlation, and spatial coherence 
and phase spectra analyses of different line parameters in the 
first paper of this series (Puschmann et al. 2003, hereafter Paper 
I), the aim of the present work is to gain insight into the varia- 
tion of physical quantities characterising granulation (like tem- 
perature (T), line of sight velocity (VlosX gas pressure (P g ) 
and density (p)), obtained from the inversion of our time se- 
ries of one-dimensional slit spectra. The present investigation 
consists of three parts. Firstly, a study of coherence and phase 
spectra analysis has been carried out, to gain information about 
the height variation of physical quantities of structures of dif- 
ferent size. Secondly, we aim at the creation and analysis of a 
model of an average granular cell, showing the granule to in- 
tergranular lane stratification of physical quantities at different 
optical depths and geometrical heights to retrieve information 
about the variation of physical quantities of different structures. 
Thirdly, the horizontal variations of physical quantities for the 
cases of a specific small and a typical large granular cell at dif- 
ferent optical depths and geometrical heights are presented. 

2. Observation 

The data used for this study consist of a 50 minute time series of 
high resolution quiet-granulation CCD spectrograms, contain- 
ing three photospheric Fe I lines (A 6494.98, 6495.74, 6496.47 
A, hereafter referred to as Line I, Line II, and Line III, respec- 
tively). Line II has been excluded from this analysis due to a 
blend of a terrestrial water vapour line. Blends in the left wing 
of Line I and in the right wing of Line III have been considered 
in the spectral synthesis, hereafter referred to as Line IV and V, 
respectively (see Table 1). 

The data were obtained on July 8, 1993, with the 70 
cm Vacuum Tower Telescope at the Observatorio del Teide 
(Tenerife). The observations have been performed at the so- 
lar disk centre. The entrance slit of the spectrograph was set 
to 100 jU. The pixel size in the focal plane of the spectrograph 
was 0"093 in the spatial direction and 2.06 mA in the spectral 
direction. In the spatial direction, the slit covered a total length 
of 95 "2 (1024pix) although after flat-fielding, only the cen- 
tral part between the two reference hairs for positioning (46 "6 
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Table 1. Atomic parameters used for the inversion of the spectroscopic time series. The columns denote in order: Line num- 
ber, element, wavelength (Fel lines: Nave et al. 1994; Ball line: Pierce & Breckinridge 1973), damping enhancement factor, 
excitation potential (Nave et al. 1994), oscillator strength (evaluated in this paper, see text for details), quantum numbers of the 
lower and upper levels of the transition (Nave et al. 1994), velocity exponent (a) and line broadening cross-section (a c ) of the 
collisional broadening kindly provided by L.R. Bellot Rubio. Line II has been excluded due to a strong blend with a water vapour 
line. 



Line 


Element 


A (A) 


E 


^exe (eV) 




Transition 


a 


<7 c (cm 2 ) 




I 


Fel 


6494.9805 


1.0 


2.39 


-1.29 


3H6.0-5G5.0 


0.247 


8.98800e- 


15 


III 


Fel 


6496.4666 


1.0 


4.77 


-0.56 


3D 2.0 - 3D 2.0 


0.279 


2.59000e- 


14 


IV 


Fel 


6494.4999 


1.0 


4.71 


-1.2 


3D 3.0 - 5P 2.0 


0.0 


0.0 




V 


Ball 


6496.9095 


0.1 


0.60 


0.1 


2D 1.5 -2P 0.5 


0.0 


0.0 





equivalent to 501 pix) was kept for further analyses. From the 
computed power spectra we estimate the effective spatial reso- 
lution achieved as ~ 0"5. 

For further details about observation and data reduction see 
Paper I. 

3. Inversion of spectroscopic data 

The inversion technique SIR (Stokes Inversion based on 
Response function, Ruiz Cobo & del Toro Iniesta 1992) pro- 
vides the stratification of thermodynamic and dynamical pa- 
rameters like temperature (T), line of sight velocity (Vlos) an d 
micro-turbulent velocity (^ m ; c ) vs. continuum optical depth at 
5000 A . Hereafter we refer our optical depths T to the standard 
wavelength A =5000 A (t = Tsooo). The geometrical height 
scale (z), electron pressure (P e ), gas pressure (P g ), and den- 
sity (p) are also evaluated under the assumption of hydrostatic 
equilibrium (HE) in each iteration but they depend strongly on 
the boundary conditions. Since a perturbation of these parame- 
ters does not significantly affect intensity profiles, the boundary 
condition can be changed afterwards. 

We refer the reader to Ruiz Cobo & del Toro Iniesta ( 1992) 
for a detailed description of the inversion method. In summary, 
the technique is a Marquardt non-linear least-square method 
(see Press et al. 1986), in which the so-called response func- 
tions of intensity are used as derivatives of the % 2 -merit func- 
tion. The procedure starts with the estimation of a model atmo- 
sphere which is modified iteratively until the synthetic profiles, 
calculated with the new model atmosphere, match the observed 
ones. The difficulty of the inversion of the radiative transfer 
equation is precisely its non linearity (in contrast to other in- 
versions like in helio-seismology, which allow the evaluation 
of the kernels once and they are kept constant throughout the 
inversion process). The response functions depend so strongly 
on the model that they have to be calculated in each iteration. 
Furthermore the iterative process has to be treated carefully 
to ensure convergence. Ruiz Cobo & del Toro Iniesta (1992) 
demonstrated the convergence of the code to reasonable and 
unique models, even in case of simulated observations with 
very poor signal to noise ratios. 

The basic assumptions of SIR are local thermo-dynamical 
equilibrium (LTE), a homogeneous plane-parallel atmosphere 
and hydrostatic equilibrium. The required absorption coeffi- 
cients, Planck function, damping and Voigt function are cal- 



culated by using Wittmann's (1974) routines modified so that 
the whole set of partial derivatives is also obtained. 

3. 1. Atomic line parameters 

The reliability of model atmospheres determined from inver- 
sions depends strongly on the accuracy of the atomic param- 
eters available for each line. For the inversion of our data set 
the values presented in Table 1 have been considered and are 
discussed below. 

Collisional broadening by neutral hydrogen atoms is known 
to be underestimated in spectral syntheses using the van der 
Waals interaction potential due to its asymptotic nature. In 
order to compensate for the smaller broadening, the micro- 
turbulence needs to be overestimated. The introduction of the 
so-called damping enhancement factor (E) in order to increase 
the broadening by collisions has the objective of mitigating this 
problem. However because of the ad-hoc nature of this correc- 
tion factor, the micro-turbulence from observed spectral lines 
is expected to be questionable. Fortunately, the use of E is no 
longer necessary following the important advances achieved on 
the quantum formulation of the collisional broadening (O'Mara 
1976; Anstee & O'Mara 1995, Barklem et al. 2000, and refer- 
ences therein), which make a realistic estimation of the inter- 
action potential with neutral perturber atoms. Therefore for the 
inversion we have used a new version of SIR (Borrero & Bellot 
Rubio 2002) implementing this new treatment of collisional 
broadening. The necessary parameters a (velocity exponent) 
and <7 C (line broadening cross-section) needed to calculate the 
collisional broadening of our lines have been kindly provided 
by L.R. Bellot Rubio (see Table 1). 

Abundances have been taken from Thevenin (1989) and the 
quantum numbers of the lower and upper levels of the transi- 
tion, central laboratory wavelengths and the excitation potential 
(fexc) from Nave et al. (1994). 

The adequate values of log gf (presented in Table 1) have 
been derived from the best-fit between observed and synthe- 
sised profiles (see Fig. 1). 

3.2. Depth grid 

For the inversion, the optical depth grid has been chosen as 
equi-spaced on a logarithmic scale, including 37 points with 
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Fig. 1. Plot of the observed and synthesised profiles using the 
atomic parameters given in Table 1. Upper panels: Mean ob- 
served profiles (dotted) and synthesised profiles (solid) and the 
difference (small panels) between them. Lower panels: Same 
but for an individual case. 



a step width of 0.2, spanning from log T= 1.2 to —6.0. The 
various nodes in which the physical quantities are evaluated 
by SIR cover the same range of optical depth. The highest 
and the deepest points of the grid are always included in the 
set of nodes. The whole iterative process is divided into sepa- 
rated groups of iterations, each characterised by the number of 
nodes selected. One of these groups runs iteratively until con- 
vergence is reached, i.e., until the variation of % between any 
two successive iterations is of no significance. The resulting 
model atmosphere enters, as an initial one, the next group of 
iterations with a larger number of nodes. By this procedure, it 
is possible to slowly decrease the number of degrees of free- 
dom. Therefore, intermediate models are smooth. For the in- 
version of our data, 3 iterations have been chosen, fixing the 
number of nodes for T to 2, 5, 7 and for Vlos to 2, 4, 6 respec- 
tively. Electron pressure P e , gas pressure P g , and geometrical 
height z are derived from the T stratification at each iteration 
assuming hydrostatic equilibrium (HE). Micro-turbulence has 
been assumed to be constant with depth. Additionally a macro- 
turbulence velocity has been derived for each pixel. 

3.3. NLTE-e ffects 

SIR calculates the atomic level populations under the LTE ap- 
proximation. Nevertheless, departures from LTE have proved 
to be important for Line I (see Paper I). A detailed investiga- 
tion of these effects has recently been presented by Shchukina 
& Trujillo Bueno (2001). The departures affect mainly the 
line source function, which deviates significantly from Planck's 
function as one proceeds outwards in the atmosphere. We have 
used departure coefficients j3] ow (t) and j3 up (t) obtained by 



Shchukina & Trujillo Bueno for the quiet sun model of Maltby 
et al. (1986 MACKKL model) and a 63 level model atom. j3i ow 
(j3 up ) stands for the ratio between the population of the lower 
(upper) atomic level evaluated in NLTE and LTE conditions, 
respectively. 

NLTE effects have to be considered for both synthesis 
of the line profiles and calculation of the response function. 
The results will be more reliable, the more similar the final 
model (obtained by the inversion) and the MACKKL model 
are. Deviations of the final model from the MACKKL model 
in deep layers do not produce significant errors, because the 
departure coefficients at those layers are always near unity for 
solar models. However, in high layers deviations of the final 
model from the MACKKL model could produce significant 
changes in the departure coefficients (see Paper I). 

3.4. Resulting stratifications of physical quantities as 
functions of optical depth and geometrical height 

The 3D box of intensity 7(x;, fj, Ak), where ; = 1, 500 with steps 
of Of'093, 1, 152 with steps of 20 sec, k= 1, 1024 with steps 
of 2.06 mA has been inverted and therefore for each spatial 
and temporal point the stratifications of the physical quanti- 
ties T, Vlos> P%< z an d p as functions of the continuum optical 
depth at 5000 A (T5000) have been obtained. The uncertainties 
of the derived physical quantities are calculated from the dif- 
ference between the observed and synthetic spectra, following 
the standard non-linear least squares error propagation. Due to 
large error bars at optical depths smaller than log 1 = — 4, only 
the optical depth range between 1 .2 and — 4 is used for fur- 
ther computations. Since the observations are obtained at disc 
centre, Vlos will be termed hereafter V z - vertical velocity. 

Subsequently the 3D boxes, containing the spatial and tem- 
poral variation of physical quantities vs. log 1 have been filtered 
layer by layer of oscillations by applying a subsonic filter in the 
Fourier space (k x , co), described in Paper I. We have chosen the 
maximum phase velocity admitted by the filter being 5 km s _1 
and a gradual cutoff of a cosine bell between the lines co = v p k x 
and co = 1 .08 v p k x has been used. 

As an example, in Fig. 2 the filtered 3D boxes of T (log T, x, t ) 
and V z (log T, x, t) are presented. Note that at log T = the struc- 
tures are almost identical to those of the intensities presented in 
Paper I. 

The geometrical scale resulting from the inversion of the 
data is different for each spatial and temporal point. The cri- 
terion to establish a common height scale has been to neglect 
horizontal acceleration at one specific height, that is, impos- 
ing horizontally constant gas pressure at this layer. A constant 
gas pressure value of Pf = \.%2x 10 5 dyncm" 3 (MACKKL 
model) at a height of zo = — 80 km (log T = 1 .08) has been cho- 
sen. Subsequently the needed shift 8z(x, t ) has been evaluated 
so that 



P % {z = z^8z)=P i 



M 



(1) 



Finally, each physical quantity T (z, x, t ), Vlos {z, x, t ), 
P g (z, x, t), p (z, x, t) has been interpolated to a new equi-spatial 
grid within a z range from — 120 km to 400 km, with a step 



K. G. Puschmann et al.: Time series of high resolution photospheric spectra 



/^^^ T (log T,X,t) 



V (log T.x.t) 



% 1 



Fig. 2. Resulting 3D boxes of the physical quantities T and V z 
as functions of logTsooo, x , t after filtering oscillations. T and 
V z are represented, with log X running from to — 4 (step width 
of 0.2), x running from 1 to 500 (step width of 0."093) and t 
running from 1 to 152 (step width of 20 sec). 



width of 20 km. The geometrical level z = corresponds to 
T5000 = 1 ■ 

4. Spatial coherence and phase spectra analysis 

An analysis by means of spatial coherence and phase difference 
spectra (see Paper I for a detailed description of this method) 
between temperature (T) and vertical velocity (V z ) fluctuations 
at the consecutive levels throughout the solar photosphere has 
been used to study the height variation of physical quantities of 
structures of different size. 

4. 1. Physical quantities in optical depth 

Coherence and phase shifts vs. spatial frequency u be- 
tween temperature and vertical velocity (i.e., between 
the pairs [T(logT = 0)-V z (logT = i)] for i = 0, -1 , — 2, -3, 
[V z (log x = -0.4) - V z (log t = /)] for i = -0.8 , - 1 , -2, -3, and 
[T (log t = ) -T (log T = z')] for i = -0.4, -1,-2,-3) resulting 
from temporal average, are presented in Fig. 3 for unfiltered 
and filtered data, respectively. 

To facilitate the estimation of the different penetration 
heights of T and V z depending on the structure size, the max- 
imum spatial frequency (minimal size) of structures contribut- 
ing to the temperature and velocity fluctuations at each opti- 
cal depth are presented in Fig. 4. This figure has been con- 
structed by assuming the following constraint: We consider 
only spatial frequencies u<2 arcsec -1 (structures > 0"5), a 
coherence > 0.5, and a phase difference of 180° ± 20° between 
temperature and vertical velocity at the different optical depths 
i, whereas between the velocities a 0° ± 20° phase difference is 
assumed. Between temperatures a 0° ± 20° phase difference at 
layers below log t = — 1 (thin dotted line of Fig. 4) and a phase 
difference of 180° ±20° at layers above (thick dotted line of 
Fig. 4) is considered. As the lower and upper limit of granular 
scales we assume ~ 0"5, as the expected limit of the spatial 
resolution of our data set, and ~ 4", respectively. 

For filtered data in general a larger penetration height for 
all physical quantities is found compared with unfiltered data. 
From the coherence and phase shift analysis between the pair 
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Fig. 3. Phase difference and coherence vs. spatial fre- 
quency for unfiltered (thin lines) and filtered data (thick 
lines) between the pairs: [T(logT = 0)-V z (logT = ;)], 
where i = 0,-1, -2, -3; [V z (logT=-0.4)-V z (logT = ;)] 
and [r(logT = 0)-r(logT=/)], where z' = -0.8,-l,-2,-3 
and i = — 0.4, — 1 , — 2, — 3 , respectively. 



[T (logT = 0)-r (log T= /)] at lower layers a 0° phase differ- 
ence together with a coherence > 0.5 at granular scales is de- 
tected (see Fig. 3 and thin dotted line of Fig. 4). A fast decay of 
temperature fluctuations with decreasing optical depth together 
with a decreasing contribution of small structures is revealed, 
vanishing at log T = — 1 . At higher layers and at granular scales 
> 1 "5 an inversion of the temperature fluctuations is observed, 
reflected in the phase jump from 0° to 180° together with a co- 
herence > 0.5 (see Fig. 3 and thick dotted line of Fig. 4). These 
findings confirm the results obtained after an analysis of global 
correlations, spatial coherence and phase spectra between the 
fluctuations of different line parameters presented in Paper I. 

For vertical velocities, the coherence and phase 
shifts between the pairs [T (logT = 0)-V z (logr = /)] and 
[V z (log T = 40) - V z (log T = i)] reveal a similar size distribution 
of velocity and temperature structures with height, as expected 
from the overshoot scenario. Vertical velocities penetrate into 
the highest layers under study, even at logT= —2.8 convective 
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Penetration of granulation — logr (unfiltered) Penetration of granulation-logr (filtered) 




Fig. 4. Maximum (minimum) spatial frequencies (sizes) of structures contributing to the T and V z fluctuations at different optical 
depths, as obtained from spatial coherence and phase difference analyses vs. spatial frequency between the pairs [ T (log T = 0) — 
V z (logT = i)L [V z (logT=-0.4)-Vz(logT = /)] and [r(logT = 0)-r(logT=/)], i = -0.2, . . . , -4, A/ = 0.2, for unfiltered (panel 
a) and filtered data (panel b), respectively. See the text for a detailed description. 



velocities of the larger granular structures (~ 4") are present. A 
less efficient height penetration of smaller structures is found, 
although still at log T ~ — 2 structures > 1 '.'A are detected. 

4.2. Physical quantities in geometrical height 

Properties observed at equal optical depths (e.g. logr = 0) do 
not refer to a unique geometric depth and observers should be 
cautious in their interpretations about temperatures, Doppler 
velocities, etc. This is a result of the rapid increase in the 
H opacity with temperature, so one looks down to shallower 
depths in hotter regions and does not see the high-temperature 
gas, because it is very opaque. As an example, much smaller 
temperature fluctuations are observed at equal optical depths 
than at equal geometrical levels (see Figs. 7 and 9). Keeping 
this in mind, the presentation and analysis of physical quanti- 
ties at different geometrical heights facilitates the interpretation 
of the results. 

Coherence and phase shifts vs. spatial frequency be- 
tween temperature and vertical velocity (i.e. between 
the pairs [T(z = km) -V z (z = /)], [V z (z = 40 km) - V z (z = i)l 
[T (z = km) - T (z = i)], resulting from temporal averages, are 
presented in Fig. 5, for unfiltered and filtered data, respectively. 
i denotes the consecutive layers in geometrical height through- 
out the solar photosphere. 

The maximum spatial frequency (minimal size) of struc- 
tures contributing to the temperature and velocity fluctuations 
at each geometrical height level are presented in Fig. 6. This 
figure has been constructed like Fig. 4, described in the previ- 
ous section. Between temperatures a 0° ± 20° phase difference 
at layers below 180 km (see Fig. 6, thin dotted line) and a phase 
difference of 180° ± 20° at layers above (see Fig. 6, thick dot- 
ted line) has been considered. As for optical depth, as the lower 
and upper limit of granular scales we have assumed ~ 0"5 (as 
the expected limit of the spatial resolution of our data set) and 
~ 4", respectively. 

Again for filtered data in general a larger penetration height 
for all physical quantities is found compared with unfiltered 
data. From the coherence and phase shift analysis between 
the pairs [T (z = 0km)- T (z = i)] the fast decay of tempera- 



ture fluctuations together with a decreasing contribution of 
small structures with height is revealed. At the layers below 
150 km and at granular scales, we find a coherence larger than 
0.5 together with a phase shift of 180° between temperature 
and vertical velocity, whereas temperatures at the different lev- 
els are in phase. At layers higher than 200 km again, as seen 
in optical depth, an inversion of temperature fluctuations for 
structures > 1 '.'5 can be observed, but only in case of filtered 
data, reflected in a phase jump between the pair [T (z = km) - 
T(z = i)] from0° to 180°. 

For vertical velocities, the coherence and phase shifts be- 
tween the pairs [T (z = 0km)- V z (z= /)] and [V z (z = 40km)- 
Vz(z = i)] confirm the results found for optical depth. The 
largest granular structures (~4") are still observed at a height 
of 370 km, whereas structures larger than 1 '.'A3 contribute until 
a height of 280 km. 

Espagnet et al. (1995) performed a phase and coherence 
spectral analysis between intensity and velocity fluctuations 
as measured at several positions (wavelengths) throughout the 
profile of the NaD2 line (X 5890) which correspond to various 
heights in the photosphere ranging from the continuum level 
up to ~ 550 km. In contrast to our results they conclude that 
temperature fluctuations of granulation do not penetrate higher 
than about 60 - 90 km. By comparing the intensity fluctuations 
at the continuum level with those at different levels in the pho- 
tosphere, they find high coherence for large granules ( > 1 !'4) 
up to 60 km, vanishing at 90 km. For Nesis et al. (1988) and 
Komm et al. (1990) the coherence of granules > 1"4 vanishes 
at 170 km above the continuum level. 

There are discrepancies in the literature about the small- 
est size of velocity structures reaching the 200 km level above 
the surface. The values range from 0."9 (Deubner 1988) to 1"5 
(Nesis et al. 1988) and up to 2"of Durrant & Nesis (1982). At 
this height, we obtain significant coherence for structures larger 
than ~ 1"3. Salucci et al. (1994) and Espagnet et al. (1995) 
claim that the velocity fluctuations remain coherent across the 
whole photosphere. Espagnet et al. (1995) compared the veloc- 
ities at different layers with both the continuum intensity and 
the velocity at the height of 30 km. From these coherence and 
phase shift analyses they conclude that the overshooting veloc- 
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Fig. 6. Maximum (minimum) spatial frequencies (sizes) of structures contributing to the T and V z fluctuations at different ge- 
ometrical heights, as obtained from spatial coherence and phase difference analyses vs. spatial frequency between the pairs 
[T{z = 0km)-V z {z = i)], [V z (z = 40 km) - V z (z = j)] and [T (z = 0km)-T {z = i)], i = 20, ... ,400km, A/ = 20 km, for unfiltered 
(panel a) and filtered data (panel b), respectively. See the text for a detailed description. 



ities of granules > 1 '.'A and > 1 '.'6, respectively, cross the whole 
photosphere up to a level of at least 550 km. As stated above, 
we detect only structures <~ 4"contributing to the velocity field 
observed at z ~ 370 km. 

The present work, in contrast to the abovementioned stud- 
ies found in the literature, is based on the application of an 
inversion method to estimate the penetration height of granu- 
lation. This allows a precise determination of the height where 
the observed perturbation of physical parameters takes place. 
Traditional methods, based on contribution functions or on the 
calculation of bisectors, can introduce large errors in the height 
determination (Sanchez Almeida et al. 1996 and references 
therein). 



5. Model of an average granular cell 

Taking into account all pixels of our time series, we have 
76 000 points available, where for each of them its model at- 
mosphere is known, i.e. the stratification of T, V z , P g , p vs. 
z as well as vs. log T. In order to obtain a model of an aver- 
age granular cell these model atmospheres have been grouped 
into 76 bins, where each bin is the average of 1000 individ- 
ual model atmospheres. As a criterion for this binning we have 
chosen the temperature values of each bin at the continuum 
layer log T = and z = 0, respectively. Thus, the resulting mod- 
els of our average granular cells, both in log T and z, consist 
of 76 averaged model atmospheres (bins), where the centre of 
the granule and its surrounding intergranular lane have been as- 
signed to the hottest and coolest bin, respectively, whereas bins 
with temperature values in between correspond to the gradual 
transition granule - intergranular lane. This is an assumption 
which might not be always fulfilled (i.e. exploding granules), 
but it seems reasonable to apply this hypothesis for the case of 
an average of granular cells. Evidently our model is not repre- 
sentative of the case of exploding granules. However exploding 
granules represent only a small part of granules; according to 
Mehltretter (1978) their number density is around 4% and ac- 
cording to Namba ( 1986) exploding granules cover about 2.5% 
of the observed area. 



5. 1. Results as a function of optical depth 

Figure 7 shows the resulting granule to intergranular lane T - 
T, V z -%, (P g -P~ g )IP~ g and (p-p)/p stratification vs. logT. 
The presented error bars are calculated by the inversion code 
for each physical quantity and individual model using the re- 
sponse functions (for a recent explanation and application of 
this concept see Socas-Navarro 2004). The error was prop- 
agated by assuming a Gaussian distribution until the values 
presented in panel e) to h). The error is minimal where the 
response functions (i.e. the sensitivities) are maximal for the 
spectral lines under study. The presentation of the physical 
quantities in panel a) to d) has been restricted to those layers 
with a reasonable degree of confidence. 

In optical depth, the absolute temperature fluctuations 
A T (granule - intergranular lane) of the average granular cell 
(panel a) are rapidly decreasing from a maximum A T ~ 400 K 
(at log T = 0.2) towards a minimum near logT= — 1, and then 
increase again in the higher layers. At about log T= — 1.5 the 
granule appears cooler than the intergranular lane, thus con- 
firming an inversion of the temperature contrast in higher lay- 
ers, which is in agreement with our results of global correla- 
tion and spatial coherence and phase shift analyses presented 
in Paper I and with the results presented in Section 4. It con- 
firms other results reported in the literature (e.g. Deubner 1988; 
Collados et al. 1996; Rodriguez Hidalgo et al. 1996, 1999) and 
it has been predicted by theoretical models (Steffen et al. 1989; 
Stein & Nordlund 1989, 1998; Gadun et al. 1997, 1999). 

The absolute velocity fluctuations (AV Z ), presented in Fig. 
7 (panel b), reach 1 .4 km s~ 1 at log T = — 0.4 decreasing slowly 
with height and are ~ at log T = — 3. Near the surface layers, 
velocities of the coolest intergranular lane are higher than the 
hottest granule ones. 

The relative pressure difference between the granule and 
the intergranular lane (AP g /P g ) decreases rapidly with decreas- 
ing optical depth and becomes practically at log T ~ — 1 (see 
Fig. 7, panel c). The hydrostatic equilibrium equation on the 
optical depth scale reads: 
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Fig. 5. Phase difference and coherence vs. spatial frequency for 
unfiltered (thin lines) and filtered data (thick lines) between the 
pairs: [T (z = 0km)-V z (z= i)] and [V z (z = 40km)-V z (z = /)], 
where i= 80, ... ,400 km, A/ = 80 km; [T (z = 0km)- T (z = /)], 
where i = 40, . . . , 200 km, Ai = 40 km. 



where g stands for gravity and k tot stands for total continuum 
absorption coefficient per mass unit. We see that the gas pres- 
sure depends only on the absorption coefficient per mass unit. 
In the Sun k tot is dominated by H opacity which, under the 
HE constraint, grows with T in lower layers and slightly di- 
minishes at higher ones (see Fig. 8). The derivative of k tot vs. 
T has been calculated for the average model. To evaluate the 
derivative, we introduced a small r -perturbation and recalcu- 
lated the ^-stratification and the absorption coefficients, im- 
posing hydrostatic equilibrium (HE). Since granules are hotter 
at deeper layers the corresponding k tot is larger and following 
Eq. 2 at larger optical depths the gas pressure in granules be- 
comes smaller compared with intergranular lanes. At the higher 
layers (beyond log T ~ — 1) small positive values of AP g /P g are 
observed as a consequence of the changing sign of both the T- 
fiuctuation and the derivative of k tot (see Fig. 7 and Fig. 8). 

The relative density fluctuations diminish faster than those 
of temperature but slower than the pressure ones. Considering 
the centre of a granule and calculating the scale height of A T/T 
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Fig. 7. Model of an average granular cell: T—T (panel a), V z - 
V z (panel b), (P g -P g )/P g (panel c) and (p-p)/p (panel d) 
stratification vs. logTsooo- Sequence from the centre of a gran- 
ule (black colour, negative velocities) to the central part of its 
surrounding intergranular lane (red colour, positive velocities). 
Colours in between show the gradual granule-intergranular 
lane transition. Panel e) to h): Corresponding error bars pre- 
sented for the centre of the granule and the central part of its 
surrounding intergranular lane. 



(hereafter H T ), AP g /P g (H P ), and Ap/p (H p ) at the deepest 
layers, i.e. the interval in log T where the relative differences 
of these parameters is lower by a factor of e, we obtain values 
of H T = 0.76, Hp = 0.48, and H p = 0.54. The intermediate be- 
haviour of density can be explained in terms of the equation of 
state taking into account that the mean molecular weight (jj.) 
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Fig. 8. Derivative of the total continuum absorption coefficient 
(k to t) vs. T, imposing hydrostatic equilibrium (HE). 



practically does not change at all. Hence at the layers where 
APg/Pg is vanishing Ap/p ~ —AT/T. 

5.2. Results as a function of geometrical height 

In Fig. 9 the resulting granule to intergranular lane T -T, V z - 
Vz> {Pg-Pg) I P% an d (p-p)/p stratifications vs. z are repre- 
sented. The error bars have been calculated by a Monte-Carlo 
simulation, e.g. by adding random noise to each physical quan- 
tity of the individual models in log 1 with a a equal to the esti- 
mated error. For each model, perturbed in this way, the cor- 
responding stratification in z is evaluated. By repeating this 
process in a large number of simulations, the resulting error is 
given by the deviation with respect to the mean of the obtained 
results. 

In geometrical height, the maximum AT increases signifi- 
cantly up to HOOK at z = - 40km (compared with the 400 K 
at log T = 0.2 in optical depth). This is in agreement with the 
results found by Rodriguez Hidalgo et al. (1999). Also accord- 
ing to numerical simulations by Stein & Nordlund (1998) there 
should be a much wider spread in temperatures at a given ge- 
ometric height (near the surface) than there is on a local opti- 
cal depth scale. Near the surface the energy transport switches 
from convective below the surface to radiative above the sur- 
face. The temperature gradient on an optical depth scale around 
logT = corresponds to an extremely steep gradient on a geo- 
metric height scale because of the extreme temperature sen- 
sitivity of the dominant H opacity, so that a small increase 
in temperature produces a large increase in opacity and hence 
a large increase in optical depth over a very small geometric 
height range. However, we obtain in geometrical height still a 
much smaller T contrast compared with observed continuum 
intensity maps since our analyses are based on a time series of 
one dimensional slit spectra and thus we are not able to cor- 
rect the data for the influence of degradation by stray light con- 
tamination. We find the T fluctuations decrease rapidly with 
height and reaching a minimum at ~ 170 km, close to the value 
obtained by Rodriguez Hidalgo et al. (1999), who localised a 
vanishing of AT at 150 km. In agreement with the results ob- 
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Fig. 9. Model of an average granular cell: T - T (panel a), V z - 
Vz (panel b), (P g -Pg) l~P % (panel c), (p -p) /p (panel d) strat- 
ification vs. z. Sequence from the centre of a granule (black 
colour, negative velocities) to the central part of its surround- 
ing intergranular lane (red colour, positive velocities). Colours 
in between show the gradual granule-intergranular lane tran- 
sition. Panel e) to h): Corresponding error bars presented for 
the centre of the granule and the central part of its surrounding 
intergranular lane. 



tained from the coherence and phase shift analysis (Fig. 5), a 
weak inversion of the T contrast at higher layers is found. 

Concerning vertical velocities, the maximum AV Z amounts 
1.4kms _1 at z = 60km (close to the value found for optical 
depth at log T = —0.4). AV Z decreases slowly from there on with 
increasing height, amounting to ^0.4kms _1 at z = 400km. 
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Fig. 10. Model of an average granular cell. Panel a): logf- 
logT stratification vs. z. Sequence from the centre of a gran- 
ule (black colour, negative velocities) to the central part of its 
surrounding intergranular lane (red colour, positive velocities). 
Colours in between show the gradual transition from a granule 
to its surrounding intergranular lane. Panel b): Corresponding 
error bars presented for the centre of the granule and the central 
part of its surrounding intergranular lane. 



This is also in agreement with results by Rodriguez Hidalgo 
et al. (1999), who found a value of ~ 1 kms -1 at z = 150 km 
and very slowly decreasing afterwards. Initially hotter and less 
dense up-flowing plasma continues to rise even beyond the den- 
sity inversion at z ~ 30 km, which clearly indicates observa- 
tional evidence for substantial overshoot into the photosphere. 
As observed already in optical depth, the velocities of the 
coolest intergranular lane are higher than the hottest granule 
ones near surface layers, in concordance with numerical sim- 
ulations by Stein & Nordlund (1998). Also, panel f) depicts 
clearly that the lines are insensitive to velocities at the deep- 
est layers. The error of vertical velocity increases dramatically 
for z < 20 km, thus the representation and analysis of velocity 
has been restricted to layers beyond. Nevertheless, by looking 
at the reliable layers, the granular and intergranular velocities 
show extrema (in absolute values) at z ~ 80 km and z ~ 20 km, 
the one of granular velocities shifted towards higher layers. 
This is in qualitatively good agreement with the simulations 
by Stein & Nordlund (1998, Fig. 5). 

The AP g values obtained (Fig. 9) could be strongly influ- 
enced by the criterion we have used to find a common geo- 
metrical height grid, where arbitrarily a layer of constant gas 
pressure was imposed. The exponential decrease of P g and the 
fact that the layer of equal P 2 has been selected deep in the 
photosphere (at z = — 80 km) can force spurious pressure differ- 
ences from local hydrostatic equilibrium at some heights in our 
model. This problem will be addressed in a forthcoming paper 
by creating a dynamic model imposing the equations of conti- 
nuity and motion. However, we can assume that temperatures 
and velocities are only weakly affected by this problem, show- 
ing a much smaller height dependence than pressure fluctua- 
tions. Densities again could be strongly affected by the changes 
introduced to correctly solve the equations of motion. The cor- 
responding panels of Fig. 9 have to be treated carefully. 

The decrease of the optical depth range (hereafter A log T) 
with geometrical height implies that granules are much less 
transparent than intergranular lanes in layers below 50 km and 
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Fig. 11. T(t ,x,log T = 0) image (panel a), corresponding T bi- 
nary image (panel b), and T binary image showing the selected 
areas (marked by rectangles) for the study of the variations of 
physical parameters across a specific small and a specific large 
granular cell (panel c). 

only slightly less transparent in layers above, see Fig. 10 where 
the log T — log T stratification is presented. Greater opacity for 
granules is observed in all layers which can be interpreted as 
follows. At layers where the granular temperature is lower than 
the average, the derivative of the continuum absorption coeffi- 
cient is also negative (see Fig. 8). 

We suggest that this change of opacity enhances somewhat 
the crossing of temperatures (inversion of temperature contrast) 
between granules and intergranular lanes around log T = — 1 
(see Fig. 7). 

6. Detailed analysis of the spatial variations of 
physical quantities across a small and a large 
granular cell 

For a comparison with the more general results obtained by 
means of studies of the average granular cell and the spatial 
coherence and phase spectra analysis, the spatial variations of 
physical quantities, i.e. for temperature and vertical velocity, 
across a small (size < 1"5) and a large (size ~ 3") granular cell 
(see Fig. 1 1 for the selected areas) at different optical depths 
(geometrical heights) are presented in Fig. 12 as temporal av- 
erages. 

Granular cells have been identified by means of a binary 
map of temperature at the logT = level, which has been 
created by applying a criterion similar to the intensity crite- 
rion for the selection of granular and intergranular regions, 
described in Paper 1. Granules (intergranular lanes) are 
defined by all those pixels in space (x) and time (f) with tem- 
perature values at log T = 0, satisfying the following conditions: 

Granules: 

r(x,f,lo g T = 0)>r(x,r,logT = 0) (3) 
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Intergranular lanes: 

r(jc,/,logT = 0)<T(jc,/,logT = 0) (4) 
where T is the mean temperature value over the whole image. 

We are dealing with time slices, thus it is difficult to really 
identify granular cells in the 'images', i.e. a small bright point 
in the x-t image could be also a part of a larger granule moving 
in and out of the slit. Likewise, it is difficult to cut through 
the centre of a large granule. The results presented here serve 
only for a qualitative comparison with the results obtained for 
the average cell. Analyses based on real 2D spectrometric data 
will provide more detailed insight in the future. By keeping 
the limitations in mind the main results can be summarised as 
follows: 

The high temperature sensitivity of the H opacity pro- 
duces much smaller temperature fluctuations at a certain level 
in optical depth than in geometrical height (seen also in the case 
of the average granular cell in Figs. 7 and 9). E.g., for the large 
cell in deep layers a AT of ~ 1800K is observed, whereas in 
optical depth only a AT of ~ 500 K is detected. 

In general, we find the variation of physical quantities much 
more pronounced across the large granular cell, which is also 
in agreement with the results obtained in Paper I by means of 
a detailed analysis of the variation of line parameters in space 
and time for a specific small and a specific large granular cell. 
In the case of the small cell at deep layers a maximum AT of 
~ 200 K in optical depth and of ~ 800 K in geometrical height 
is found. For the larger cell a very broad T maximum at the 
centre of the granule is found, also in agreement with the results 
of Paper I. AT decreases rapidly from its maximum of 500 K 
with decreasing optical depth (i.e. from 1 800 K with increasing 
geometric height). The decrease occurs even faster in the small 
cell. 

For vertical velocities, as in the case of the average model 
(see Figs. 7 and 9), we find similar values in optical depth 
and geometrical height. For the small cell, the maximum AV Z 
amounts to ~ 1.4kms , whereas for the larger cell a maxi- 
mum AV Z of ~ 2 km s _1 is found. 

In contrast to the average case (see Fig. 9) no larger veloc- 
ities are found at the intergranular lane as compared with the 
granule. Note that in Fig. 9 the largest downflows appear only 
for the coolest bin, which do not coincide with the samples pre- 
sented here. 

We also point out an asymmetric distribution of velocities 
in particular cases, e.g. at deep layers of the larger cell, the 
maximum up-flow is located near the granular border whereas 
at higher layers and in the case of the smaller cell the vertical 
speeds show increased blue shifts towards the granule centre. 

This is in agreement with the results obtained in Paper I 
by comparing the line core velocities of the two Fe I lines A 
6494.98 A and A 6496.47 A. We conclude that in particu- 
lar cases of very large granules an asymmetric distribution of 
vertical velocities across the cell can be observed. This is in 
agreement with the findings of other authors, e.g. de Boer et al. 
(1992), Nesis et al. (1993), Hanslmeieret al. (1994), Hirzberger 
(2002), indicating that maxima of up-flows in large granules do 



not always coincide with the maxima of intensity, and velocity 
peaks adjacent to the intergranular lanes can be observed. 

7. Conclusions 

The vertical structure of the solar photosphere has been inves- 
tigated by means of different approaches. In Paper I, an analy- 
sis of global correlations, spatial coherence and phase spectra 
between different line parameters, obtained from a time series 
of ID spectrograms, has been carried out. In the present work 
the spatial and temporal distribution of the thermodynamic al 
quantities and the vertical flow velocity is derived as a func- 
tion of optical depth (logr) and geometrical height (z), apply- 
ing the inversion technique SIR (Stokes Inversion based on 
Response functions) on the above data set. Spatial coherence 
and phase spectra between the fluctuations of temperatures and 
vertical velocities vs. optical depth and geometrical height have 
been computed for unfiltered and filtered data to determine the 
height variation of physical quantities of structures with differ- 
ent size. A model of an average granular cell has been derived, 
showing the granule-intergranular lane stratification of differ- 
ent parameters as a function of log T and z- The cases of a spe- 
cific small and a specific large granular cell have been investi- 
gated to compare with the results obtained in the more general 
cases. 

Regarding the granular temperature pattern, a fast decay of 
the temperature fluctuations between granules and intergranu- 
lar lanes (AT) with increasing height can be observed. A less 
efficient penetration of smaller cells to higher layers is reflected 
and already at ~ log T=— lorz~170 km the T - T coherence 
is lost at all granular scales and Ar shows a minimum. At the 
layers beyond, for structures > 1 '.' 5 an inversion of the temper- 
ature contrast in both optical depth and geometrical height is 
found. We show how the effect of opacity somewhat enlarges, 
in the logr scale, the inversion of the temperature contrast, 
which has been ascribed so far only to a cooling produced by 
adiabatic expansion. Furthermore, in agreement with the nu- 
merical simulations by Stein & Nordlund (1998), the high tem- 
perature sensitivity of the H opacity produces at equal optical 
depths much smaller temperature fluctuations (we report a Ar 
of 400 K at log T = 0.2) than at equal geometrical levels (where 
we find a Ar of 1400 K at z = — 40 km). A maximum Ar of 
200 K (800 K) in optical depth (geometrical height) for a spe- 
cific small granular cell and of about 500 K (1800K) for the 
larger one, decaying rapidly with height, confirm the results 
obtained for the average cell. 

A much slower decay with height of the vertical convec- 
tive velocity fluctuations between granules and intergranular 
lanes (AV Z ) is found. In deep photospheric layers, a maximum 
AV Z of 1.4kms _1 in both optical depth and geometrical height 
is in concordance with the 1.2kms~' and 2kms~ 1 found for 
the specific small and specific large granular cell. Going up 
in height, as in the case of temperature, a less efficient pen- 
etration of convective velocities of smaller cells is reflected, 
though even at log T ~ — 2.0 and z = 280 km structures > 1 "4 
are found. At z^ 370 km, only velocities of structures at the 
largest granular scales (~ 4") are present. This similar size dis- 
tribution between velocities and temperatures with height to- 
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Fig. 12. Panel a) to d): Spatial studies of the variation of T and V z across a particular small and large granular cell at different 
optical depths. Panel e) to h): The same for geometrical heights. Different colours (sequence from black to red) are indicating 
different depths (heights). T is presented from log T = 0.2 (z= — 40km) up to logT= — 3 (z = 400km) , V z from log T= — 0.2 
(z = 20 km) up to log X = - 2.6 (z = 360 km and z = 400 km). 



gether with the temperature inversion clearly provides observa- 
tional evidence for substantial overshoot into the photosphere, 
far beyond the density inversion, which we find near the surface 
layers. At deep photospheric layers, the behaviour of the verti- 
cal velocities reflected in the simulations of Stein & Nordlund 
(1998) is for the first time qualitatively reproduced by obser- 
vations. In our model the centre of the granule shows smaller 
velocities than the centre part of the intergranular lane. Both 
velocities (in absolute values) reach extrema near the surface 
layers, but the granular one is shifted towards higher layers. 
Concerning the specific studies of a small and large granular 
cell we point out an asymmetric distribution of vertical veloci- 
ties for the large granular cell at the deep layers, where the max- 
imum up-flow is located near the granular borders. However, at 
higher layers and for the specific small cell the vertical speed 
shows a more symmetric distribution with increased blue shift 
towards the granule centre. This confirms findings of other au- 
thors, e.g. Hirzberger (2002), indicating that maxima of up- 
fiows in large granules do not always coincide with the max- 
ima of intensity, and velocity peaks adjacent to the intergranu- 
lar lanes can be observed. 

The obtained gas pressure differences between granules 
and intergranular lanes (AP g ) could be influenced by the cri- 
terion we have used to find a common geometrical height grid, 
defining arbitrarily a layer with horizontally equal P 2 . A deep 
photospheric layer has been selected for the horizontal pressure 
equilibrium. This together with the exponential decrease of P 2 
can introduce spurious pressure differences from hydrostatic 
equilibrium at some heights in our model. This problem will 
be addressed in the next paper of this series, by imposing the 



equations of continuity and motion to our model, thus creating 
a 2D semi-empiric dynamic model of an average granular cell. 
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